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ABSTRACT NOMENCLATRUE 

A transient two-phase numerical model is T Temperature (K) 
developed to investigate insulation performance of t Time (s) 
cryogenic propellant tanks using finite difference method. Cp Specific heat at constant pressure (J/KgK) 
The model includes the mathematical formulations for aes Mass transfer across interface (Kg/s) 

. oe d K Thermal conductivity (W/m k) 
convective and radiative heat transfer from the ambient, 

P Tank Pressure (Pa) 


heat and mass transfer across liquid-vapor using two- 


; : ; : Gr — Grashoff Number 

phase formulations. The model is validated with the Pr Prandil Number 
experimental test data on LH2_ tank. Subsequently, Re Reynold Number 
analyses are carried out to determine heat in-leak, g Gravity (m/s2) 
evaporation rate, evolution of pressure and temperature in h Coefficient of heat transfer (W/m2K) 
ullage of cryogenic tank with different foam insulation Q Heat transfer rate (W) 
thicknesses of 20mm, 25mm, 30mm and 35mm. Analysis q Heat flux (W/m2) 
shows that tank pressure rise is significantly higher at A Surface area of heat transfer (m2) 
lower insulation thickness. The present study brings out a H Enthalpy (/kg) 

E : D Tank diameter (m) 
correlation between rate of tank pressure rise and foam ; 
; P ; : : X Ullage height (m) 
insulation thickness which shows that LN2 evaporation : ; 

aa i xX Insulation thickness (mm) 

rate has negligible effect on tank pressure rise for <i velocity (m/s) 
simulated conditions. A minimum foam insulation V Volume (m°) 
thickness required for no ice formation over foam outer 
surface is also suggested in the study. Greek Symbols 
Keywords: Foam Insulation, Heat In-Leak, Cryogenic, HM Dynamic viscosity (Pa-s) r say 
Tank Pressure, Liquid Nitrogen o Stefan-Boltzmann constant (5.67*10°W/m‘K’) 


Coefficient of Thermal Expansion (K") 
Absorptivity 

Density (Kg/m’) 

Parameter for tank geometry (m’') 
Parameter for thermal properties (Pa-m’/s) 


e<evewnw 


Subscripts 
ull Ullage 
int Interface 
liq Liquid 
Ww Wall 
ins Insulation 
f Fluid 
a Ambient gas 
(e) Outer insulation surface 
I Inner insulation surface 


INTRODUCTION 


The proper design of insulation system on a cryogenic 
propellant tank is vital for the nominal operation of 
cryogenic engine. The significant external heat load would 
affect the propellant thermal conditions during the mission 
which would affect the engine performance. Flight 
cryogenic tank uses foam insulation. Proper thickness of 
foam insulation is required to maintain the propellant 
temperature within allowable limits and also avoids tank 
pressure rise beyond a limit during the entire mission. 
Cryogenic engine requires minimum tank pressure at the 
start which can also be met by properly designing 
insulation thickness. Insufficient insulation can lead to 
possible ice formation on outer tank surface, which would 
increase the lift-off mass resulting payload penalty. Thus, a 
detailed thermal analysis of cryogenic tanks insulation 
system is essential to achieve nominal performance of the 
mission. 

Several investigations were reported in literature for 
insulation of cryogenic tanks. Wexler [1] proposed a 
model which describes the evaporation rate of liquid 
helium in cryogenic tank. A heat transfer study in 
cryogenic tank was presented by Khemis et al [2]. The 
numerical model developed by Khemis et al [2] showed 
large disparity with their experimental data. Later a brief 
review was given by Khemis et al [3] through the 
experimental study. Subsequently, Boukeffa et al [4] have 
experimentally evaluated the heat losses in cryogenic 
cryostat. 

A numerical heat transfer study considering 
convective flux over the cryogenic tank was carried out by 
Pictet [5]. Arsonval [6] proposed the insulation system for 
tank using double wall with vacuum annular section. This 
tank was later improved by Dewar [7] to store hydrogen 
liquid at 20.4 K. Several experimental studies were carried 
out by the researchers [8-12] to investigate the heat 
transfer in cryostat from ambient. Shu et al.[13] 


recommended for an optimal insulation of a nitrogen 
cryostat with a number of screen layers. 

Though several numerical studies are carried out to 
determine heat transfer in cryogenic tank with insulation, a 
comprehensive numerical model for determining the effect 
of foam insulation system on tank thermodynamics are 
unavailable. Thus, the present study aims at bringing out 
the impact of foam insulation thickness on heat in-leak, 
tank pressure, fluid temperature and ice formation over the 
tank surface. The model uses a finite difference approach 
to predict thermo-fluid transients at liquid-vapor interface 
of cryogenic propellant tank. 


FORMULATION AND SOLUTION METHODOLOGY 


A transient two phase flow and thermal model is 
developed in SINDA/FLUINT simulator. The model is 
used to analyze the evolution of temperature and pressure 
inside the tank considering conduction, convection, 
boiling, condensation, and heat and mass transfer across 
the interface. A schematic of the propellant tank showing 
various heat transfer mechanism is depicted in Fig.1. 





























FIGURE.1 SCHEMATIC OF PROPELLANT TANK WITH 
DIFFERENT HEAT TRANSFER MODES 


A cylindrical tank of 0.5m diameter and 1m height is 
taken for the analysis. The tank is foam insulated and its 
thickness varies as 20mm, 25mm, 30mm and 35mm which 
is the investigation parameter for the analysis. The tank is 
80% filled with saturated LN2 at lbar pressure. Outer 
boundary of the foam is provided with solar and ambient 
heat fluxes. Transient analyses are carried out for 15min 
duration, within which the tank pressure rises due to 
external heat loads. The model uses mathematical 
formulations for heat transfer from ambient to tank surface 
by convection and radiation and heat transfer from the 
ullage gas to the wall as well as that to the liquid nitrogen. 
Three-dimensional conduction, with — temperature 
dependent material properties, is modeled for the tank 
walls and foam insulation layers. Thermodynamic state 
equation is used to solve for fluid conditions. Iterative 
analysis is carried out with various insulation foam 
thickness of 20mm, 25mm, 30mm and 35mm. The effect 
of insulation thickness on various parameters like tank 


pressure, ullage temperature, heat in-leak and evaporation 
rate are investigated which are the key parameters used to 
ensure the thermal requirements for an optimum insulation 
thickness. 

The heat transfer modes considered in the analysis are 
expressed mathematically as below: 


Heat transfer from ullage to liquid 


Heat, from ullage to bulk liquid, transfers through 
interface. However, the heat transfer from ullage to 
interface and interface to liquid are equal. 


Qin-tig = Q.n-int i Qini-tig (1) 


Heat transfer from the interface to liquid is expressed 
as: 


On-tig = PogetiqAint Tint = Tig) (2) 


The heat transfer from the ullage gas to interface is 
due to convection. Since, diffusing or condensing 
constituent carries its own individual enthalpy, there is an 
additional heat transfer across interface due to mass 
transfer of these species. Effective heat transfer from the 
ullage gas to interface is expressed as: 


Qynint = Men-intAine Tir — Tine) + MA 3) 


m and H are mass and enthalpy of nitrogen transfer 
across interface. m is positive for condensation and 
negative for vaporization of the liquid and can be 
expressed as 


[Pine tg Aint Tine ~ Tig) — 
7 H 


It has been assumed that the heat transfer is due to natural 
convection with the heat transfer coefficient being 
expressed by the empirical correlation [14]: 


int-lig/*int int ull-int u 


Aint (T; i Ti) 
(4) 





k, 
h=K,C—X : (5) 
l 
Where: 
X =(Gr)(Pr) (6) 
and 
a 
Hy 
Pr=| 2A (8) 
K, 


Constant (C) is 0.54 and 0.27 for heat transfer from 
ullage to interface and interface to liquid respectively. n = 


0.25 and Ky (heat transfer adjustment factor) is set to 1.0. 
The length scale (1,) is set to the diameter of the tank. 


Heat transfer from wall to Ullage 


The heat transfer from wall to ullage gas is expressed 
as: 


Qa = Avy (T,, = T 1) (9) 


It has also been assumed that the heat transfer is due 
to natural convection and heat transfer coefficient is 
expressed as: 


k 
h =0.508 7 Pr°* (0.952 + Pr) °* Gr (10) 


Ss 


ky is thermal conductivity of fluid and 1, is height of 
tank at particular point where h is calculated. 


Heat transfer from ambient 


Convection and radiation models are used to calculate 
incident heat on the tank. Coefficient of convective heat 
transfer on tank surface is calculated by using Churchill 
and Bernstein [15] formulation based on wind velocity of 
2m/s. Heat transfer through radiation is calculated by using 
solar heat flux data (solar) and altitude dependent ambient 
temperature value. Total heat transfer from outside is 
express as: 








On: = Os = Qua d 1) 
Ooi = Aeon Peay (T er a De ) (12) 
0.62Re °° Pr®33 sig 4/5 
ho =|0.3+— 1+{ Ree ¥o (13) 
conv [1+(0.4/Pr)*?)?” 28200 D 
V,-D 
Rea (14) 
MH, 
Qvad = GiiAaty a COA an 5 i. =3 spa ) (15) 


To calculate heat-in-leak from ambient to the tank 
wall, three dimensional transient heat conduction equations 
are solved considering the variation of thermal 
conductivity of the insulation with temperature. Three- 
dimensional transient heat conduction equation is given by 
following expression considering K as a function of 
temperature: 


OT OT OT _ Cp OF 
ax? dy? @z?—s«K (TT) Gt 


All the above mentioned heat transfer modes 


(16) 





represent the thermal mathematical model which have 
been incorporated in the analyses. 


RESULTS AND DISCUSSION 


Transient two-phase thermal analysis is carried out to 
determine the effectiveness of foam insulation for various 
thermodynamic parameters such as tank pressure, ullage 
temperature, evaporation rate and heat in-leak. The model 
is first validated using LH2 tank pressure data from launch 
vehicle. Subsequently, the mathematical model is applied 
to study the effect of foam thickness on tank 
thermodynamic parameters and foam surface icing. Results 
of the analysis are discussed below. 


Mathematical Model Validation 


The numerical model is validated with the 
experimental data obtained from the test on LH2 tank. 
Corresponding tank dimensions are considered for 
validation. Figure 2 shows the comparison of the 
numerical model with the experimental test data of LH2 
tank pressure during ground hold condition. 
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FIGURE.2. MODEL VALIDATION WITH 
EXPERIMENT DATA 


Steady pressure rise of 0.7 bar is measured in 60s 
during test. It can be seen in Fig.2. that the pressure rise 
from the numerical model is well corroborated with the 
measured test data. The mathematical model is further 
used for insulation effectiveness studies as discussed 
below. 


Effect of insulation thickness on tank pressure 


Tank pressure evolution: Tank pressure 
evolution during the 900s transient analysis for various 
foam thicknesses are shown in Fig.3. 
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FIGURE 3.TANK PRESSURE WITH DIFFERENT 
INSULATION THICKNESS 


It can be seen from the figure that pressure rise is 
maximum for 20mm insulated tank and reduces with 
thicker tank insulation. This is because the tank pressure 
depends largely on ullage temperature which increases due 
to external heat sources. Poor insulation allows more heat 
into the tank leading to higher pressure rise. Pressure rise 
beyond tank structural limit is fatal for cryogenic storage 
tanks. Hence, better insulation is essential for safe 
operation of cryogenic tanks. 


Rate of pressure rise: Figure.3 shows maximum tank 
pressure rise for 20mm insulation thickness. From the 
figure, it is also evident that the rate of pressure rise is 
different for each insulation thicknesses. Since rate of 
pressure rise is steady close to 900s, a study is carried out 
to find a correlation between the rate of tank pressure rise 
and insulation thickness. 


Pressure rise in cryogenic tank is due to temperature 
rise as well as liquid evaporation. LN2 evaporation has 
two effects on tank pressure; raising it through ullage mass 
addition and dropping it by cooling the ullage. In the 
following mathematical formulation, pressure rise due to 
evaporation is assumed negligible and ambient heat 
sources is assumed to affect the pressure rise significantly, 
which gets convected into the ullage through conductive 
heat transfer across insulation thickness. Pressure rise rate 
in ullage based on temperature rise rate can be expressed 
as: 





(dP) 41 = OR (AT) 41 (17) 


dt dt 


The heat balance across the foam can be written as: 


(AT) ins (I, -T;) 





Ca a mCp a f KP mai (18) 
G7) 
KinsAy ul xX = My un pean ~ Tu) (19) 
(dT), 
hy Ay—un (T,, = Ten) a MitCP ut dt : (20) 


where, X is the insulation thickness. Since heat 
transfer path for ullage temperature rise is conductive 
through insulation, rate of pressure rise can be 
parameterized from Eqns. 17, 18, 19 and 20 as follows: 


(AP) _ (Bee . Dif Ay ut } i 








dt Cp Vin J X 
(dP), 1 
——" =(¢).(w).— (22) 
it (p).(Y) = 
Where 
k. R(T -T, 
p= ins ( oO i) (23) 
Cp 
y — Ay uit (24) 
Vin 


@ is thermal dependent parameter assuming steady 
pressure rise with (T,-T;) constant, whereas yw is tank 
dimension dependent parameter. 


y= Z + = For cylindrical ullage 
x D 
(D/x) 





For spherical ullage 


~ [(D/2)—(x/3)] 


Where, x is tank ullage height and D is the tank 
diameter. 


Thus, pressure rise parametric equation independent 
of tank dimension can be written from Eqn. 22 as 


dP 1 1 
( dull —= [¢].— (25) 
dt yw xX 


From the Eqn 25, it is evident that pressure rise rate 
follows inverse power trend with insulation thickness. For 
the simulation conditions, y value is found to be 173.56. 


Subsequently, the pressure rise parameter in Eqn 25 is 
analysed for different insulation thicknesses of 20mm, 
25mm, 30mm and 35mm, considering LN2 evaporation in 
the model, and results are shown in Fig. 4. 
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FIGURE 4.TANK PRESSURE RISE RATE WITH 
DIFFERENT INSULATION THICKNESS 


The analysis results match within +3% of the 
correlation for rate of pressure rise given in Eq. 25, which 
was formulated assuming negligible LN2 evaporation rate. 
Variation in mathematical correlation and simulation 
results could be attributed to evaporation been considered 
in simulation. This shows that effect of rate of evaporation 
of LN2 on tank pressure rise for the simulation conditions 
is negligible and pressure rise in the tank is mainly due to 
heating of ullage by external sources. 


The inverse power correlation stated in Eq. 25 shows 
that higher insulation thickness has lesser effect in rate of 
pressure rise whereas reducing the thickness of insulation, 
significantly increases the rate of pressure rise, which is 
undesirable. 


Effect of insulation thickness on evaporation rate 


Different foam thicknesses are taken to study its effect 
on rate of LN2 evaporation during the transient analysis. 
Variation of LN2 mass evaporation rate, per liquid surface 
area, for different insulation thicknesses is shown in Fig.5. 
LN2 upon heated up from external sources gains latent 
heat energy which leads to evaporation. From the plot, it is 
evident that liquid evaporation rate reduces as tank 
insulation thickness increases. In cryogenic propellant 
tanks, evaporation of propellant results in reduction in 
useful propellant mass thereby, reducing engine thrust 


duration. Hence, it is vital to reduce liquid evaporation in 
cryogenic tanks by increasing the insulation thickness. 














en de 

E 

2 10 mi 

g 8 Wig 

i w 

5 6 he 
E 

5 4 

S 

uu 2 


15 20 25 30 35 40 


Insulation Thickness (mm) 


FIGURE 5. EVAPORATION RATE WITH DIFFERENT 
INSULATION THICKNESS 


Effect of insulation thickness on ullage temperature 


The study of evolution of ullage temperature for 
20mm, 25mm, 30mm and 35mm insulation thicknesses on 
the tank is carried out and results are presented in Fig.6. 
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FIGURE.6 GAS TEMPERATURE IN ULLAGE WITH 
DIFFERENT INSULATION THICKNESS 


It is to be noted that the graph is plotted from 
liquid-ullage interface to the top of the tank. The interface 
is always at saturation temperature and gradient in ullage 
temperature starts afterwards as seen in Fig.6. As the 
insulation thickness reduces, maximum value of ullage 
temperature increases as seen in Fig.6. Liquid surface 
temperature also increases which is equivalent with 
saturation temperature rise corresponding to tank pressure 
increment. 


Effect of insulation thickness on foam temperature 


The analysis is carried out to study the foam 
insulation temperature profile, at 0.1m from top of the tank 
which is at ullage side, for different insulation thicknesses. 
Variation in foam temperature pattern for different 
insulation thicknesses at 900s are shown in Fig.7. 
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FIGURE 7.FOAM TEMPERATURE ALONG THE 
THICKNESS 


It is evident from the figure that higher insulation 
thickness results in higher temperature at outer surface of 
the foam insulation. This shows that the ambient heat-in 
leak will be minimum for higher insulation thicknesses as 
convective heat transfer from ambient is directly 
proportional to temperature difference between foam outer 
surface and ambient. In another perspective, temperature 
gradient within the foam is highest for 20mm insulation 
which means conductive heat transfer through the foam is 
highest, resulting in maximum heat in-leak. 


Effect of insulation thickness on external heat in-leak 


Heat in-leak: Analysis is carried out to study 
heat in leak variation with different insulation thicknesses 
and results are shown in Fig.8. 
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FIGURE 8. HEAT-IN-LEAK INSIDE THE TANK WITH 
DIFFERENT INSULATION THICKNESS 


The results show that heat in-leak reduces with 
increase in insulation thickness. This result is expected as 
foam insulation temperature profile shown in Fig.7. 
denotes higher temperature at foam surface for higher 
insulation thickness, resulting in low ambient convective 
heat transfer, as discussed in previous section. 


Ice formation on tank surface 


Water vapor present in ambient condenses and freezes 
over the tank insulation outer surface, when temperature at 
the outer surface of the insulation drops below 273K. This 
results in deposition of ice layer over the foam surface, and 
may peel off during ascend phase of launch vehicle. Icing 
phenomenon leads to undesirable increment of lift-off 
mass and its peel-off may result in damage of surrounding 
hardware. Analysis is carried out to determine the foam 
outer surface temperature, at 0.4m from tank base which is 
at liquid side, for different foam thicknesses. Solar heat 
flux is considered in one case and not considered in other 
to analyze worst case scenario. Foam outer surface 
temperature variation with foam insulation thickness is 
shown in Fig.9. 
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FIGURE.9 FOAM OUTER SURFACE TEMPERATURE 
WITH DIFFERENT INSULATION 
THICKNESS 


It is evident from the figure that icing will not be seen 
when solar heat flux is incident on the tank. But, in case 
the tank is place in shade where little or no solar heat flux 
is incident, possibility of ice formation over foam surface 
increases as foam insulation thickness reduces. Lower the 
insulation thickness, lower will be the foam outer surface 
temperature which would raise the ambient heat in-leak. 
Figure 9 shows that for thickness of insulation less than 


30mm, foam surface temperature is lower than 273K and 
thus, ice will be formed over the surface. Hence, insulation 
thickness of 30mm or above is required to avoid icing over 
foam insulation. 


CONCLUSION 


Transient two-phase thermodynamic model of 
cryogenic LN2 tank is developed to investigate the effect 
of foam insulation thickness on tank parameters such as 
pressure, ullage temperature, heat in-leak and liquid 
evaporation rate. Validation of mathematical model is 
carried out with experimental test data on LH? tank. 
Subsequent transient analysis shows that increasing the 
insulation thickness reduces the heat in-leak and thereby 
decreases tank pressure rise which is crucial for flight 
cryogenic tanks. Ullage temperature and liquid evaporation 
rate are found to be lower when higher insulation 
thicknesses are used. Pressure rise rate for different 
insulation thicknesses are mathematically analysed to 
bring out a correlation independent of tank geometry. The 
mathematical correlation when compared with analysis 
results shows that rate of LN2 evaporation has negligible 
effect in rate of pressure rise for the simulation conditions 
considered. The correlation also shows that for lower 
insulation thickness, rate of pressure rise is significantly 
higher. Effect of presence and absence of solar heat flux 
incident on the tank is also determined for different 
insulation thicknesses. When little or no solar heat flux is 
incident on the tank, foam outer surface temperature drops 
below 273K for insulation thicknesses less than 30mm, 
leading to formation of ice over the foam surface. To avoid 
its adverse effect, insulation thicknesses of 30mm or above 
are recommended. 
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